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ABSTRACT 


Spherically symmetric perturbations in AdS-scalar field systems of small amplitude 
£ approximately periodic on time scales of order 1/e 2 (in the sense that no significant 
transfer of energy between the AdS normal modes occurs) have played an important role 
in considerations of AdS stability. They are seen as anchors of stability islands where 
collapse of small perturbations to black holes does not occur. (This collapse, if it happens, 
typically develops on time scales of the order 1/e 2 .) We construct an analytic treatment 
of the frequency spectra of such quasiperiodic perturbations, paying special attention to 
the large frequency asymptotics. For the case of a self-interacting 0 4 scalar field in a 
non-dynamical AdS background, we arrive at a fairly complete analytic picture involving 
quasiperiodic spectra with an exponential suppression modulated by a power law at large 
mode numbers. For the case of dynamical gravity, the structure of the large frequency 
asymptotics is more complicated. We give analytic explanations for the general qualitative 
features of quasiperiodic solutions localized around a single mode, in close parallel to our 
discussion of the probe scalar field, and find numerical evidence for logarithmic modulations 
in the gravitational quasiperiodic spectra existing on top of the formulas previously reported 
in the literature. 



1 Introduction 


Numerical investigations in the field of non-linear stability of Anti-de Sitter (AdS) space- 
times [1,2] have revealed a patchwork of globally regular and collapsing solutions [3-5]. In 
the corresponding non-linear perturbation theory [6-10], one finds a matching competition 
between pro- and anti-stability features [7—9] . In particular, there are powerful sets of 
selection rules [7,8,11,12] limiting the efficiency of energy transfer between the AdS normal 
modes in the weakly non-linear regime, both for scalar fields coupled to gravity, as in [ -10], 
and probe self-interacting scalar fields in non-dynamical AdS, as in [13,14]. 

Ultraviolet dynamics is of crucial importance to the problem of evolution of AdS per¬ 
turbations. Instability and black hole formation occur through non-linear energy transfer 
to very short wave-length modes, thus the question of ultraviolet asymptotics in weakly 
non-linear perturbation theory comes to the fore; see, for instance, the recent paper [10]. 

An important aspect of the non-linear perturbative AdS dynamics is the existence of 
periodic and approximately periodic solutions. For example, an algorithm to construct 
exactly periodic solutions order-by-order in the field amplitude has been presented in [4] 
(and extended to massive fields in [15]). More broadly, one can consider quasiperiodic 
solutions, in which the AdS normal modes oscillate with constant amplitudes at leading 
order in the field amplitude on time scales inversely proportional to the square of the field 
amplitude [ ]. Since these time scales are precisely the ones on which collapse occurs in 
numerical experiments, as well as the shortest time scales on which substantial energy 
drift may occur in the weakly non-linear regime in AdS settings [6,7], the existence of such 
solutions presents an obstruction to efficient collapse. Islands of stability in the phase space 
of AdS perturbations are thought to be anchored by such quasiperiodic solutions [9,16]. 

At leading order, effectively resummed perturbation theory [6-8] takes the form of equa¬ 
tions describing slow drifts (on time scales inversely proportional to the square of the 
field amplitude) of the amplitudes and phases of the AdS normal modes. This leading 
order approximation can be derived by many techniques (multi-scale or renormalization 
group analysis, averaging) and goes by a few different names in the literature (‘two-time 
framework,’ ‘renormalization flow equations,’ ‘time-averaging,’ ‘resonant approximation’). 
Quasiperiodic perturbations appear in this context [6] as solutions for which the normal 
mode amplitudes remain constant in time, despite the presence of weak non-linear inter¬ 
actions. This condition can be reduced to an infinite system of cubic algebraic equations, 
solutions of which were analyzed numerically in [6], revealing exponential damping for high 
frequency modes. It was pointed out in [ ] that the abundance of quasiperiodic solutions 
is intimately related to selection rules in perturbation theory, which, in turn, give rise to 
conserved quantities [8]. The solutions were explicitly parametrized by conserved quantities 
in [9] and numerically shown to be stable in [16]. 

Here, we commence our studies of the ultraviolet behavior of the AdS perturbation the¬ 
ory by analyzing the ultraviolet asymptotics of the spectra of such quasiperiodic solutions. 
We do it both for self-interacting probe scalar field dynamics [ ] and for gravitational 
interactions [6]. For the self-interacting probe scalar field, we develop a fairly complete 
analytical treatment that recovers the essential features of the asymptotic spectra in a way 
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that agrees with direct numerical evaluations. Specifically, one gets amplitude spectra ap¬ 
proaching an exponential decay multiplied by the mode number for high-frequency modes. 
For the gravitational case, both the numerical behavior and the structure of the equations 
are more involved. We display the asymptotic properties of the coefficients in the equa¬ 
tion determining the quasiperiodic spectrum, and explain the general qualitative features 
of the solutions, such as the existence of two spectral branches previously pointed out nu¬ 
merically [6]. We are unable, however, to draw clear analytic conclusions on the spectral 
asymptotics from these equations. At the same time, our numerical study of the gravita¬ 
tional quasiperiodic spectra gives evidence for a hierarchy of slow logarithmic modulations 
superposed on the simple fitting formulas mentioned in [6]. 

We present our results below by first reviewing (in section 2) the improved perturbation 
theory in AdS at first non-trivial order, leading to an accurate description of the dynamics 
on time scales inversely proportional to the square of the field amplitude, as well as the 
construction of quasiperiodic solutions in this context. We then describe our analysis of 
the quasiperiodic spectra asymptotics, first for a self-interacting scalar field (section 3), 
and then for the gravitational case (section 4). Other aspects of ultraviolet asymptotics 
of AdS dynamics are discussed in our parallel work [17], where, in particular, we analyze 
the ultraviolet behavior of the interaction coefficients of the time-averaged equations in 
arbitrary dimension. 

2 Time-averaged effective theory and quasiperiodic solutions 

We shall briefly review the construction of effective flow equations describing the slow 
energy transfer between the AdS scalar field modes, starting with the simpler set-up of a 
self-interacting scalar field [13, I ] in AdS^+i. A considerably more systematic discussion 
of these matters (phrased in a context including full gravitational non-linearities) can be 
found in [7,8] . 

We start with the global AdS^+i metric written as 

ds 2 = ——— (— dt 2 + dx 2 + sin 2 x dfl 2 ,,) , (1) 

cos 2 x v ' 

where the radial x coordinate runs between 0 and 7t/2 and dn 2 l _ 1 is a (d — l)-sphere line 
element. The equations of motion for a spherically symmetric probe scalar field t) with 
a 0 4 /4! non-linearity take the form 


- d 2 (j) + -— lz]—d x (tan d 1 x d x <t>) = 
tctn x 


3! cos 2 x 

Expanding this in normal modes 4>(x,t) = J2 n c n (t)e n (x) yields 
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C n + UJ n C n | 
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where 


e n (x) = k n (cosx) d P r 


d_-t d> 
9 -*■’ 9 ; 


(cos(2x)) with 


k n — 


2 yjn\(n + d — 1)! 


\n 


+ 1 ) 


( 4 ) 


P n are Jacobi polynomials, u> n = d + 2 n is the peculiar integer (fully resonant) AdS^+i 
frequency spectrum, and 


(-'ijkl. 



(x)e j (x)e k (x)ei(x) 


tan d 1 x , 

-5 —dx 

cos^ x 


( 5 ) 


are the ‘interaction coefficients’ that will play a key role in our discussions below. 

One can now zoom in on the small amplitude regime, assuming amplitudes of order e, 
and switch to the ‘interaction picture,’ introducing slowly varying complex amplitudes a k 
(which would be exactly constant if the non-linearities were neglected): 


c n =e ( a n e UjJnl + a n e 
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and hence 
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By the standard lore of averaging [8,13], when e is small, only non-oscillating terms on 
the right-hand side can contribute significantly. This implies that, out of the sum over jkl , 
only terms satisfying the resonance condition u) n = ±c Oj ± uj k ± Ui survive (with all the 
three plus-minus signs a priori independent). Furthermore, the C coefficients turn out to 
vanish (see footnote 3 of [8] for a simple proof) for all such resonance conditions other than 
c o n + ujj — ito k + u>i, i.e., n + j = k + l (and similar relations obtained by interchanging j 
with k or /). One thus arrives at the following effective time-averaged theory accurately 
approximating the original equations on time scales of order 1/e 2 : 


‘2ioj n d n —— ] (■ ’vJk.n ■ k^n - j k: 


( 8 ) 


jk 


where 3/3! = 1/2 accounts for two other contributions interchanging j with one of the two 
other summation indices in (7). 

One may look, following [6,14], for solutions of the time-averaged equation preserving 
the normal mode amplitudes |afc|, i.e., solutions in which there is no energy flow between the 
modes, despite the presence of non-linearities. 1 A large class of such solutions is given by 

1 The structure of the energy flow can be elucidated by writing (8) in terms of real amplitudes A n and 

phases B n , a n = A n e lBn , as we do for the gravitational case in (42-43). In the equation for A n , the right- 
hand side of (8) is converted to a sum of terms proportional to AjA k A n+ j_ k sin(£? n + Bj — B k — B n+ j- k ). 
Each such terms describes the influence of modes j, k, n + j — k on the amplitude (energy) of mode n via 
non-linear interactions. The condition we impose below on the phases, B n = Bq +n(B i — Bq), forces these 
terms to vanish individually, signifying the absence of resonant energy transfer within any given quartet of 
modes. 
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imposing a n = a n e * e2/3nt , where /3 n = /3 0 + n(Pi — /%), and a n , /5 0 , Pi are time-independent. 
This ansatz converts the differential equations (8) to a system of algebraic equations for 

a n ,P o,Pi- 


2ca n ci n (/3o T ti(Pi Po)) 




( 9 ) 
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It is our main objective in this paper to examine what we can say analytically about 
solutions to the quasiperiodicity conditions of this sort for the time-averaged system. For 
the self-interacting scalar held, we shall present a rather complete picture of the ultraviolet 
properties of solutions with spectra localized around a single mode. 

We shall also present some considerations for the case with full gravitational interactions 
and a free massless scalar held [ ]. The ideology of our treatment is very similar though both 
technical aspects and phenomenology are more complicated and the amount of progress we 
can make is considerably more modest. In this case, one considers a spherically symmetric 
scalar held coupled to a dynamical spherically symmetric metric of the form 


ds 2 


L 2 

COS 2 X 


i~A - 25 S ^ n2 x I 


( 10 ) 


with A and 5 being functions of x and t. Due to spherical symmetry the metric is de¬ 
termined on each constant time slice by the constraint equations in terms of the matter 
distribution (i.e., the scalar held prohle). One hence essentially obtains an effective theory 
for (j), albeit with non-linearity much more complicated than in (2). Improved perturbative 
expansions can be constructed for this system [6-8] , though doing this analytically [7, 8] 
requires considerable effort. See, in particular, [ ] for a formulation in the language of time¬ 
averaging. Once the effective averaged description of the slow energy transfer has been 
constructed, one can look for quasiperiodic solutions along the lines we described above for 
the probe scalar held. 


3 Non-linear probe scalar field in AdS 

3.1 Quasiperiodic solutions dominated by a single mode 

The conditions (9) defining the spectra of quasiperiodic solutions form an infinite system 
of cubic algebraic equations, clearly beyond direct reach of any analytic methods. 

A particular approach to constructing solutions of (9) has been proposed (in the gravi¬ 
tational context) in [6]. Namely one looks for solutions with one mode being much larger 
than any other mode. (There are obvious exact solutions with only one mode excited, and 
singlc-mode-dominated solutions can be seen as perturbations of such obvious single-mode 
solutions.) In [6], such solutions were constructed numerically. Here, we shall focus on 
solutions that are strongly localized in mode number, for which, as it turns out, a fairly 
complete analytic treatment can be given. 

Consider hrst for simplicity a solution dominated by mode 0, namely do = 1 (there is 
a scaling symmetry in the equation that always allows to set the amplitude of any given 
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mode to 1), a± = 8 <C 1. One can then see that there is a consistent ansatz a n = q n S n , 
which results in the leading powers on the right-hand-side and left-hand-side matching in 
all of the equations (9). The procedure of constructing the solutions then goes as follows. 
In each equation (9) one only keeps the leading powers of 8, matching its coefficient on the 
two sides of the equation. For equation number 0, this defines /3q , for equation number 1, 
this defines /?i, each subsequent equation determines < 72 , <? 3 , etc. Importantly, restricting 
to the leading power of 8 (strong localization limit) in each equation truncates the infinite 
sums in (9) to a finite number of terms. Specifically, 

1 x k ' 

2ukqk(/3o + k(/3i — (3 0 )) = - Ckoj,k-jQjQk-j■ (H) 

Z 3 =0 

(One of the summation indices in (9) had to be set to 0 to ensure a leading, rather than sub- 
leading, power of 5.) Equation (11) presents a tremendous simplification over (9) and can 
be analyzed quite thoroughly. We shall do that after examining the large index properties 
of the C coefficients. 

For solutions dominated by mode number J > 0 , a slighly more complicated, but still 
completely transparent picture emerges. Namely, one assumes aj — 1, aj+ 1 = <5 < 1. The 
self-consistent ansatz to substitute into (9) is then taken to be a*, = qk^ k ~ J \. With this 
ansatz, equation number J determines /3 0 + J(/3i — /3 0 ). Equations number J + 1 and J — 1 
thereafter determine /3 0 , /3i and qj-\. These equation are effectively quadratic and have two 
solutions, which explains the two branches seen in [6,11]. Proceeding further iteratively, 
each pair of equations number J + n and J — n with n < J forms a 2 x 2 linear system 
determining qj+ n and q.j- n ■ Finally, each equation with number greater than 2 J determines 
the q with the corresponding number, much like what (11) does for the J = 0 case. We 
shall see below how this procedure works explicitly for J = 1. 

We briefly comment on the general nature of the small 6 expansion, of which we are 
mainly interested in the leading terms. (Systematic investigations of this issue are out¬ 
side the scope of our present study.) The pattern we have displayed above suggests an 
asymptotic expansion around 6 = 0 of the form 

«„ = (?„ + tfs + q^S 2 + •••). ( 12 ) 

Indeed, after the leading order contribution has been determined in the way we have de¬ 
scribed, subleading terms can be fixed order-by-order by substituting (12) into (9) and 
solving the resulting linear equations. Importantly, the prefactor 8 completely drops 
out when substituted in the equations (the leading order is 5-independent, as seen, for 
example, in (11) for J = 0). For that reason, the fact that the asymptotic expansions for 
different a n start with different powers of 8 by no means upsets the consistency of deter¬ 
mining subleading corrections simultaneously for all a n . Note that, for corrections of each 
given order in 8, the infinite sums in (9) will truncate to a finite number of terms. Indeed, 
in (9) written for mode number n, the right-hand side contains an infinite sum over j and 
k, with terms bearing factors of 8^~ J \+\ k - J \+\ n +3~ k - J \_ Since all terms in the power of 8 are 
positive and the Erst two terms grow as j or k are increased for j, k > J, it is evident that, 
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for each given n and J, only a finite number of choices of j and k exist that keep the power 
of 5 below any given value. This leaves a finite number of terms out of the sum in (9) at 
each n and at each order in 5. 

Asymptotic expansions of the form (12) strongly suggest that the leading order results 
we deal with here become as accurate as one needs for sufficiently small 6. The corrections 
in (12) are uniquely determined by solving linear equations and simply slightly shift the 
existing leading order spectral branches we have described above. 


3.2 Asymptotics of the interaction coefficients 


We shall now analyze the asymptotic behavior of the C coefficients appearing in (11) for 
AdS/i. The mode functions for global AdS^+i are given by (4). For large mode numbers, 
we can use the approximation (see, e.g., [1 ]): 


in x) a+ 2 (cos xY + 2 P^’^ (cos(2a;)) 


CM - 

n 


(13) 


to write the mode functions as 

2 A/rXnWl)T(nTd) 


e n [x = 


r(n + 0 
2 y/rXrTTTjTXnTd) 

V^T (n + § + \) 


(cos x) d Pn 2 1,a ^ (cos(2a;)) 


^r^-(cotx) 2 |cos ((d + 2n)x — ^(rf — 1)j + O ^ | . (14) 


For generic dimension d, the approximate expression for the mode functions will diverge 
near x — 0, and we will not be able to employ the straightforward approach for extracting 
the asymptotics of C coefficients that we utilize below for d = 3. At d = 3, 


e n (x) « 




(15) 


does not diverge. 

We pause to note that this formula is closely related to, and could also be easily derived 
from, the exact simplified expressions for the AdS 4 mode functions that appeared in [16] 
while we were preparing the current article for publication. Those expressions were used 
in [16] to obtain closed-form (albeit very long) expressions for the interaction coeffients 
in the gravitational case, and this could also be done for the case of the probe scalar. 
While the availability of such closed-form expressions for the interaction coefficients in 
this particular dimension is certainly a welcome development, the simpler analysis of the 
asymptotics provided here will be sufficient for our purposes in the present paper. In 
passing, we also note that an efficient recursive way to compute interaction coefficients in 
arbitrary dimension is provided in [17]. 
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We shall start by displaying the asymptotic formula for C with four large indices, which 
is an instructive aside: 


C'nrnhl — 


16 f n72 sin(a; n a;) sin(o; m a;) sin^a;) sin(u;;a;) 


da; 


71 


Sill X 


— — I — UJ n + UJ m + UJk + Ui\ + \U) n — UJ m + UJk + LOi\ + CJ n + UJ m — UJk + UJi\ 

71 V 

+ \u n + + Wfc — CJ; | — \uj n + UJ m + UJk + CJi \ — \u) n + UJ m — UJk ~ UJi \ 

|^Jn ^Jm ~b ^Jk ^Jl | Y^n ^/c T til/• ( 16 ) 

In a general number of dimensions, we expect the first powers of the combinations of 

frequencies in this expression to generalize to power d — 2, in analogy to the considerations 

of [!7]. 

The actual coefficients appearing in (11) have at least one small index, while the other 
indices may become large asymptotically. For three large indices, we have 2 


4a/2 3 

eo(a;) = —— (cosa:) 


7T 


(17) 


and 


a 


nmkO —" 


32^/2 r /2 sin(c 1 ; ri a;) sin(o; m a;) sin^a^cosa;)" 


71 


32 v / 2 


x 


71 


In particular, 


da; 

sm a; 

"0 if max{n, k,l} > (n + k + l + 2)/2 

7t/ 16 if max{n, k, 1} — (n + k + l + 2)/2 

37t/ 16 if max{n, k, 1} = (n + k + l + l)/2 

^ 7t/ 4 otherwise 


o /o 

Cokj,k-j -~ 3.6013. 

7T 


(18) 


We proceed with two small and two large indices, using 

ei(x) = ^-^(cosa:) 3 (3 — 8(sina;) 2 ). 


3^ 


(19) 


( 20 ) 


a 


128a/ 3 r 72 


nmOl —" 


37T 2 


sin(a; ri a;) sin(o; m a;)(cosx) 4 (3 — 8(sina;) 2 )da; 



V/32 

if 

m — n 

128a/3 

7tt/64 

if 

77i — n 

3tt 2 '' ' 

5tt/32 

if 

771 — 71 


lo 

otherwise 


( 21 ) 


2 Here and below, we are using the explicit formulas PQ U ’ b \x) = 1 and P-) a ’^(a;) = |(2(a + 1) + (a + b + 

2)(x-l)). 
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Figure 1: Co, n + m , n ,m plotted as a function of n and m. 


So in particular C 0} i >n , n +i —> ~ 3.6755. 


128 /‘ 7r /^ 

Cnmoo——v- / sin(o; n a;) sin(u; m a:)(cosa;) 4 da; 


7T 



" 7t/64 

if 

m — n 

128 

—— x ^ 

7t/16 

if 

77i — n 

7T 2 

3vr/32 

if 

771 — 77 


lo 

otherwise 


( 22 ) 


So in particular Co,o,n,n —* ^ ~ 3.8197. 

These observations for large mode number limits taken in different directions in the mode 
number space reveal a very simple plateau-like topography, see Fig. 1. Almost everywhere 
in the large index value region, with one index kept equal to 0 according to (11), the plot 
of C is completely flat at a fixed height of 8a/2/7t. If one goes to the border regions of 
this plateau, where one of the three large indices becomes small, the plot rises somewhat, 
but neither is the rise numerically significant (only a few percent), nor is it statistically 
significant when appearing inside the sums, as in (11), since this rise is only attained for 
special, non-generic values of the indices. In the end, the large mode-number asymptotics 
of the solutions to (11) can be characterized almost entirely in terms of the single ‘plateau’ 
value of C, i.e., 8 \f2fi7. We shall see below how this works explicitly. 
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3.3 Solutions dominated by mode 0: asymptotic estimates 


Armed with asymptotic expressions for the C coefficients, we can start onr analysis of the 
large mode-number limit of the solutions to (11). As we already explained, one assumes in 
(9) ao = 1, a\ = 8 -C 1 and a n = q n 8 n . From the equations number 0 and 1 in (9), one 
obtains 

a _ CoOOO _ Cooil 

" ~4cT’ “ 2(d + 2)' 



For d = 3, this becomes 




(24) 


At higher orders in 8 , one has to solve (11). Note that the C coefficients inside the sum in 
( 11 ) are nearly constant at large k , except for the endpoints of the sum, where they vary 
slightly. Since there are considerably more terms in the interior of the sum than near the 
endpoints, it is natural to assume that the equation can be approximated well by simply 
taking C = 8 v^ 2 / 7 r (the plateau value we have identified above), solving the equation and 
then checking that our assumption has been justified. It is likewise natural to assume that 
the term proportional to k 2 in the coefficient of q k dominates the lower powers of k. All 
these assumptions can, of course, be explicitly checked by an a posteriori comparison to 
numerical evaluations of the solutions to ( 11 ). 

Applying the above simplifications to (11), one obtains: 


4\/2 i k? i 

4fc 2 (/A - f3 0 )q k = - V 'qm-j, or —q k = (25) 

7T z ' IU z ' 

3 = 1 i =1 

One can assume the simple asymptotic behavior q n = Qn 1 e Kn , and check for what values 
of the parameters k, 7 , Q it is compatible with (25). We find that hi completely drops out 
from (25). It therefore cannot be determined by simple asymptotic matching and requires 
solving (11) iteratively mode-by-mode. After hi has dropped out, one is left with 


k -1 

F +2 = I0V2 Q E f(k-jy. (26) 

3 = 1 

For large k, we can approximate the sum by an integral over ( = j/k: 


k -1 


E 


P(k- i) 1 ^ 


k W / dCC(l-C) 7 . 


(27) 


Matching the leading powers on both sides of the equation fixes 7 = 1 . Thereafter, (26) 
fixes Q = We thus expect the asymptotic spectrum to behave as 


3a/2 

a n - n (Se ) 


( 28 ) 
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Note that e K cannot be absorbed in <5 since the normalization of 5 is fixed at the infrared 
end of the spectrum by imposing a\ = 6. We have checked numerically that (28) accurately 
approximates exact solutions to (11). In particular, we have verified that no logarithmic 
modulations are present in addition to the power-law modulation in (28); we will explain 
the method we used in section 3.5. The situation will not be similarly simple when we turn 
to the gravitationally interacting case. 

3.4 Solutions dominated by mode 0: systematic analysis 

We have explained above that natural simplifying assumptions (which can be verified a 
posteriori by direct comparison with numerics, even in the absence of rigorous error bounds) 
reduce the quasiperiodicity conditions (11) for solutions dominated by mode 0 to (25). We 
have then verified that the asymptotics (28) is analytically consistent with the simplified 
equation (25), and also provides a valid approximation to the full numerical solutions of 
( 11 ). 

At a practical level, one could rest content with (28). However, it would be better to 
have a more explicit analytic control over how this asymptotics emerges. (Note that being 
asymptotically consistent with an equation does not, strictly speaking, imply that solutions 
to this equation, with initial data specified in the infrared, approach asymptotics of this 
form.) Equation (25) provides a relatively simple iterative definition of q n via a quadratic 
recursive relation reminiscent of the one for the Catalan numbers. It can be analyzed 
quite exhaustively, with a curious mathematical structure emerging, which is what we shall 
explain next. 

To gain further analytic insights in (25), we introduce a generating function 

OO 

q(X) = 1 + e A T q n e Xn 

n=2 

(25), together with the initial conditions qo = qi = 1, is then equivalent to 

q" = e~ x + 10 y/2(q - l) 2 . (30) 

This equation has an obvious mechanical interpretation: a particle is moving in the g-space 
and A-time with an external force e _A and a potential —10y / 2(g — l) 3 /3. One is looking 
for a special solution of this equation that stops at A = +oo at the inflection point of the 
potential, i.e., q — 1. One-dimensional mechanics in a cubic potential is well-known to 
display a run-away behavior in finite time (all solutions run down the potential slope to 
q = Too at a certain finite moment of time). In the vicinity of that run-away moment, 
which we shall call hi, the e _A term in (30) is negligible and one can integrate the equation 
to obtain 

3a/2 

10 (A - hi) 2 ' 

Of course, the particular value of the run-away moment hi is sensitive to the entire evolution 
and requires solving (30) exactly, which cannot be done analytically. 




(29) 
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What have we learned? q defined by (29) is an analytic function of z = e~ x with a 
double pole singularity at z — e~ K . q n are the coefficients of the Taylor series expansion 
of this function around z — 0. It is generally understood that high order behavior of 
the Taylor series is governed by the singularities of the expanded function in the complex 
plane. Darboux’ theorem (see, e.g., [19]) gives a precise statement of this correspondence. 
In particular, a double pole at z — e~ K precisely corresponds to q n ~ ne Kn . Furthermore, 
the coefficient of (31) matches the coefficient of (28). Overall, our systematic analysis of 
(25) based on generating functions completely corroborates the simple asymptotic estimates 
of the previous section. 


3.5 Solutions dominated by mode 1 

As a proof-of-concept demonstration, we shall now describe how to analyze solutions dom¬ 
inated by mode 1. Compared to the solutions dominated by mode 0 considered above, this 
case presents the feature of having two spectral branches, originally seen for gravity in the 
numerical simulations of [6]. Extension to solutions dominated by still higher modes, on 
the other hand, is a pure technicality. 

We assume a spectrum cq = 1, <22 = 5 1 and a n = q n 8^ n ~ l \ substitute it in (9) and 

retain only the leading powers of 5 in each equation. (9) with n = 0,1, 2 gives 

2wo<?oA) = 2 ( 2 Ciooi 9 o + C1102) , (32) 

2cni/3i = -Cun, (33) 

2w 2 (2/3i — /3q) = -(Cn2o<?o + 2 C 1221 ) ■ (34) 

For d = 3, /3i = 19/187T. Solving the two remaining equations simultaneously leads to 


2 185^ 7 

« + ~mT qa + 3 


0, 


(35) 


with two solutions, corresponding to two spectral branches of the type seen numerically 
in [6]: 


-185^+^42726 347 - ^7121 

Qo ~ 264 ’ _ 2527T ’ 

-185^-^42726 347 + ^7121 

Qo ~ 264 ’ “ 2527T ' 

The remaining equations reduce to 

k—1 k 

4cUfc(/3o + k(/3 1 — Po))qk — Qo^^Cq ,k,j,k-j<lj1k—j + 2 k ; l +! - ~ I ■ 

I-1 3—1 


(36) 

(37) 


(38) 
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In the limit of large k, we found that, for d = 3, the interaction coefficients are approx¬ 
imately constant: 


c. 


8y/2 


o ,k,j,k-j 


7T 


— 8 . 11(1 C\ k,j,k+l—j ~ — Ci. 

7T 


(39) 


By evoking similar logic to what we have already used in our mode 0 derivations (with¬ 
out recapitulating all the subtleties), we assume an asymptotic ansatz q n = Qri < e Kn and 
substitute it in (38), obtaining 


4 k4.(A) + k(/3i — /3o))Ad — 


k— 1 k 

QoCoQY, - dV + CiQe^ + 1 - j) 1 ■ 

3 =1 3=1 


(40) 


Matching the leading power of k on both sides of this equation results in 7 = 1. We also 
obtain the following relation between Q and k, 


48(A - fa) 

QqCq + Cie K 


(41) 


with the specific numerical values for /3 0 , /3i, g 0 , Co, C 1 given above. 

We then solve (38) numerically and fit the results to q n = Qri y e Kn { 1 + - + %). This is 
done by fitting 5 adjacent points to the above expression, and then moving these 5 points 
to higher and higher values of n in hope to see convergence. Quick convergence does, in 
fact, follow. For the solution (36), we get 7 = 1, n ~ —0.06, Q ~ 0.0846. For the solution 
(37), we get 7 = 1, n ~ —0.648, Q ~ 1.75. Both sets of values of n and Q are in agreement 
with (41). We therefore conclude that our analytic treatment of the large mode-number 
limit provides a sound perspective. Attempts to augment the asymptotics with logarithmic 
modulations, for example, a factor (Inn) 71 , and obtaining estimates for 71 , result in a 
quick convergence of 71 to 0 as the mode number increases. This indicates the absence of 
logarithmic modulations, in accordance with our analytic results. We shall now turn to the 
case of gravitational interactions, where logarithmic modulations are, in fact, present. 


4 Gravitational interactions 

4.1 Gravitational quasiperiodic solutions 

In the introduction, we have briefly explained that, compared to the self-interactive probe 
scalar held we have considered above, the effective time-averaged description for the full 
gravitational non-linearities is conceptually quite similar, but technically much more com¬ 
plicated. The averaged equations, as derived [8], are of the form (we are now using real 
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amplitudes A n and phases B n , related to the complex amplitudes by A n e lBn = a n ): 
UJ} + i k A 


2ui dAi 
e 2 dt 


2ujiAi dBi 
e 2 dt 


- ^ SijkiAiAjAh sin (B t + B k - B { - Bj), 

i j k 


( 42 ) 


Wi +t Uj —UJk 

ij=l UJ} ^ {k,l} 

TiA'i — RuA^Ai — SijkiAiAjA k cos(£>; + B k — £/ — Bj). 

i i j k 


Ldi -\~CJj —<-^k +^Z 


(43) 


The coefficients of the terms appearing in these time-averaged equations are expressed in 
terms of certain integrals of products of the AdS mode functions (4): 


Ti = l^Xuu + |y)„, + 2 u!w[°l 0) + 2 ufwg®, 


(44) 


Ra —~ 


1 f uj 2 +U}?\ , 9 „ x . „ ( ufYiM - ufYnii 


2 \uf-uf 

, ,2 , ,2 

UJ] UJ A 


(ujfXun — (jj^Xia ,) + 2 




+ -x(Yuu + Yini) 


(Xuu - Xm) + fcffcfWg” + RS") + < 45 ) 


Jijkl . \ I d" ~t~ J ( UJjpJjU} k Xnj k U>iYnj k ) 


4 yu^ + CUj COi — UJ k Ujj — co k 

1/1 1 1 

+ 


4 y cjj + ojj Ui — oj k ojj — uj k 


(ujjLJ k LJiXij k i UJjYjikl) 


If 1 


4 \uji + (Uj uji — co k ujj — u k 

1/1 1 1 


T ] (uiUkUiXjiki cjjYijki) 


(ujjUJj id i X k jjj u> k Y ik ji). (46) 


4 \uji + ujj uji — uo k ujj — uo k 
The integrals in the expressions above are defined by 

/•tt/2 

X ijk i= / d xe' i (x)e j (x)e k (x)ei(x)(]i(x)) 2 u(x), 

J o 

rn/2 

Y ijkl = / dx e' i (x)ej(x)e' k (x)e' l (x)(n(x)) 2 ufx), 

Jo 

rtt/ 2 px 

= / dxe i (x)e j (x)n(x)u(x) / dye k (y)ei(y)fi(y), 

Jo Jo 

pn/2 px 

= j d xe' i (x)e' j (x)]i{x)u(x) d ye k (y)ei(y)n(y). 


(47a) 

(47b) 

(47c) 

(47d) 
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The construction of quasiperiodic solutions [6] of (42-43) is implemented analogously 
to the procedure we have described for (9). One first assumes B n = (/3 0 + n(/3i — /3 0 ))e 2 t, 
whereupon (42) ensures that A n are constant (this is how quasiperiodic solutions of the 
sort treated here are defined). Substituting the ansatz for B n into (43), on the other hand, 
leads to an infinite system of algebraic equations for /3 0 , /3 i, A n , analogous to (9). 

One can furthermore concentrate on solutions strongly localized around one chosen 
mode, as we did in section 3.1, with the localization controlled by an adjustible parameter 
S. The resulting structure is straightforwardly analogous to (11), with the small 5 limit 
truncating the infinite sums in quasiperiodicity condition to finite sums, and solution for f3 0 , 
/3i and (iteratively) A n becoming possible. In particular, the same explanation as in section 
3.1 for the existence of two spectral branches holds. (The branches, present for solutions 
dominated by modes other than mode 0, were originally observed numerically in [6].) 


4.2 Asymptotics of the interaction coefficients 

As for the probe scalar case, the asymptotic analysis of the quasiperiodicity condition for 
gravity involves understanding the properties of the interaction coefficients appearing in 
(42-43). Of particular importance are the S coefficients with two or three large indices 
(since these are the coefficients controlling the complicated sums coupling a large number 
of different triplets of modes). 

For d = 3 one can apply a similar sort of analysis to what we described in section 3.2 to 
extract the values of the X and Y integrals when three indices becomes large. For k > 0 
and 0 ^ j A, k, one gets 


X, 


_ .-Ts? if 0<j<k 
0if J = 0,fc 


(48) 


Xkj(k-j) o 


(75+10fc)y/2 

87T 

7y/2 


if 0 < j < k 
if j = 0, k 


(49) 


Xj(k—j)kO 


(45-10j)y / 2 

87T 

13%/2 


47T 

q \/2 


if 0 < j < k 

if j = 0 

if j — k 


(50) 


YkOj(k-j) — 


32 v / 2 


x 


7T- 


661* -l^L k + 7 ~lfj(k - j) - ff k 2 + §jk(k -j) if 0 <j<k 

if j = 0, k 


256 


245* 39* j 13* l 2 

128 32 ^ 32 


(51) 


fo kj(k-j) — 


96^ f Hi — Wsk - If j(k-j) if 0 <j<k 

if j = 0, k 


x 


256 

5n 

32 


(52) 
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Figure 2: S'o.n+m^m plotted as a function of n and m, viewed from two different angles. 

Qfi /o (M-mJ + M k + ¥Ak-j) if 0 < j < k 
Y (k-j)jk0 = — X if + %k + %k 2 if i = 0 (53) 

11 ^ 3 = k 

With these values, one can compute the corresponding S coefficients. One discovers, how¬ 
ever, that the dominant contribution, which would scale like the square of the index numbers 
according to (46), in fact exactly vanishes. This leaves S growing as the first power of the 
index number (under uniform scaling of the indices) when 3 indices become large, and when 
the resonance condition is satisfied. 

One can perform similar analysis for S with only two indices becoming large. In this 
case, cancellations do not occur and one discovers that S grows like the index number 
squared (under uniform scaling of the two large indices). 

The above considerations leave one with a topography of the plot of the S coefficients 
(see Fig. 2) much more complicated than the simple plateau we have observed asymptoti¬ 
cally for the case of the probe scalar held. For example, for Sokj,k-ji if both indices k and 
j are scaled up uniformly in proportion to A, the ^-coefficient also scales in proportion to 
A. But if j (or k — j ) is allowed to stray close to small values, the plot rises dramatically, 
attaining values proportional to k 2 . 

4.3 Solutions dominated by mode 0 

We can now describe the analytic structure of gravitational quasiperiodic solutions strongly 
dominated by mode number 0. Just as we did for the probe scalar field, one assumes a 
spectrum of the form A 0 — 1, Ai — S <C 1 and A n = q n 5 n and substitutes it in (43), 
together with B n = (/3 0 + n(/3i — /3o))e 2 t, keeping only the lowest powers of 5 in each 
equation. Equations number 0 and 1 then determine /3q and 8\. The remaining equations 
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take the form 


k -1 

{2Wfc(A) + — /?o)) + -Rofc} Qk = ~ Sk-jjokQjQk-j- (54) 

j'=i 


One cannot, however, proceed with this equation in the same way we did for the probe 
scalar held case. The reason is that the dependence of S on its indices favors the summation 
index regions close to the end-points in (54), where q itself cannot be approximated by its 
asymptotics. At this stage, we cannot see any obvious ways to proceed with our asymptotic 
analysis and extract the large mode number behavior of q, even if one has detailed control 
over the asymptotics of S and R (for example, via the exact analytic expressions of [16]). 

One might be surprised that the exact asymptotics for the gravitational case proves so 
elusive, given that simple fitting formulas, q n ~ e Kn /n in our language, have been proposed 
for solutions dominated by mode 0 in [ ]. (Note that, just like for the scalar case, e Kn drops 
out if substituted in (54) and it simply corrects 5 n to the actual asymptotic exponential 
fall-off, given that the normalization of 5 is fixed in the infrared by imposing 5 = A±.) If 
the asymptotic solution is as simple as q n ~ e Kn /n, it should seem strange that one cannot 
at least verify it using the equations and encounters all the analytic difficulties we have 
observed. Given this perception, we have decided to re-examine the issue numerically. Even 
in the absence of successful asymptotic analysis, (54) provides a straightforward iterative 
way to evaluate q n , which can be easily programmed in Mathematica. We have implemented 
this procedure up to mode number 70 in hope of getting a better sense of the asymptotic 
behavior of q n . 

We first attempted to fit our numerical findings to q n = Qn r e Kn , which is what we have 
used for the scalar case. All of our fits are performed by choosing a group of adjacent points 
with the number of points equal to the number of unknowns in the putative asymptotic 
formula, fitting exactly to this group of points, and then moving the group of points to 
higher and higher index values in search of convergence. For q n = Qri y e Kn , we observed 
very slow convergence for 7 . While not being very far from —1, 7 fails to reach a plateau 
all the way until mode number 70. 

We have then attempted the same procedure with a logarithmic modulation q n = 
Qn 7 (lnn) 7 l e K " and observed that the convergence of 7 is considerably improved, while 
71 does not appear to go anywhere close to 0. (Note that a similar procedure applied to the 
probe scalar held spectrum yields 71 rapidly converging to 0.) We finally attempted to test 
the presence of In (In n) corrections, though admittedly our mode range is not particularly 
large to allow tracing very slow modulations reliably. Nonetheless, applying the fitting 
procedure to 


q n = Qn 7 (lnn) 7 l (lnlnn) 72 e' < 


c i c i 

1 H-1--- 

n n In n 


c 2 

+ ^T + 


n 2 n 2 In n 


+ 


n 2 In 2 n 


(55) 


results in good convergence within the examined mode range. 7 rapidly reaches — 1 , within 
a percent, while 7 ! and 72 appear to be around 0.3 and 0.1, respectively. (These values 
change by only a few percent between estimating around mode 50 and around mode 70 with 
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our procedure, indicating adequate convergence. Removing the corrections decreasing as 
1/n 2 and faster from the above ansatz makes the plots less flat, though the actual estimated 
values for the powers remain close to the ones quoted above.) We thus find evidence that 
quasiperiodic spectra of AdS perturbations with gravitational non-linearities display slow 
logarithmic modulations in their ultraviolet asymptotics, making them considerably more 
challenging for both analytic and numerical treatments than the power laws previously 
mentioned in the literature. It would be interesting to analyze such logarithmic corrections 
further, especially with the explicit expressions for the interaction coefficients in AdS 4 that 
have just appeared in [ 16]. The numerical optimization allowed by these expressions makes 
it possible to go to much higher mode numbers than what one can get by brute force 
evaluation. (See also our parallel work [ ] that presents effective iterative techniques for 

evaluating the interaction coefficients in a general number of dimensions.) 
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